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A  molecular  point-dxpole  model  shat  has  been  used  successfully  to  calculate 
ferroelectric  and  electro-optic  properties  of  the  perovskite  barium  titanate  (SaTiO's). 
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la-tes  lorentz  factors  by  the  method  of  Rama  Rau  and  the  crystal’s,  internal  electric 
fields  by  the  Ewald  method  -of  lattice  sums.  The- model  further  treats  the  polariza¬ 
bilities  of  ions  as  field-dependent  and  simulates  covalency  in  bonds  by  replacing  the  : 
full  ionic  charge  by  an  ''effective''  charge-  The  Kerf  and  Pockel  electro-optic 
constants  are  determined  by  the  anisotropy  induced  in  the  refractive  index  by  an 
external  electric  field-  V 

The  calculations  done  showed  that  the  model  could  match  the  experimental  refrac¬ 
tive  index  and  dielectric  constant  ana  calculate  satisfactory  electro-optic  constants 
, in  the  cubic  (paraelectric)  phase,  but  was.  unable  to  calculate  correctly  the  natural 
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Preface 


In  the  late  1940’s,  the  study  of  ferroelectricity  received  a  new 
momentum  and  direction  with  the  discovery  of  ferroelectric  barium 
titanate  (BaTiO,)  and  identifies tion  of  its  crystal  structure.  The  new 
interest  lay  in  the  fact  that  the  tetragonal  crystal  structure  of  ferro¬ 
electric  BaTiO^  was  simple  in  comparison  to  the  structure  of  other 
ferroelectrics  and  might  lend  itself  to  en  explanation  of  its  ferro¬ 
electricity  through  a  description  of  its  physical  structure.  Lead 
titanate  (PbTiO^)  was  also  found  at  about  the  same  time  to  be  ferro¬ 
electric  at  room  temperature,  with  the  same  tetragonal  structure,  and 
many  attempts  have  since  been  made  to  calculate  the  properties  of  both 
materials.  Two  treatments  that  have  been  particularly  successful  in 
calculating  the  basic  ferroelectric  properties  are  a  thermodynamic 
treatment  and  a  point-dipole  treatment.  Very  recently.  Dr.  Donald  C. 
Wunsch,  now  at  the  Air  Force  Weapons  Laboratory,  succeeded  in  extending 

the  point-dipole  model  to  calculate  the  electro-optic  behavior  of  BaTiCL  « 

5 

Further  verification  of  the  model’s  ability  to  perform  these  calculations 
correctly  is  looked  for  here,  with  an  application  to  PbTiO^  . 

I  am  indebted  to  Dr.  Wunsch  for  sparking  my  interest  in  this  field 
and  for  having  established  the  necessary  mathematical  tools  and  guide¬ 
lines  to  use  in  working  with  PbTiO^  ,  also  for  the  discussion  and  time 
spent  in  mailing  little  problems  out  of  big  ones.  X  thank  my  Thesis 
Advisor,  LtOol  John  S.  DeWitt  for  the  'unending'  patience  extended  to 


me  during  this  venture. 


August  S.  Sanches 
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TESTING  A  POINT-DIPOLE  ELECTRO-OPTIC  MODEL 
BY  AN  APPLICATION  TO  LEAD  TITANATE 

I.  Introduction 

Background 

In  a  doctoral  dissertation  entitled  "Phenomenological  Theory  of 
Refractive  Index  and  Electro-Optical  Phenomena  in  Crystals  of  Perovskite 
Structure  (With  Application  to  Barium  Titanate)",  Wunsch  (Ref  27)  has 
proposed  a  point-dipole  model  for  calculating  electro-optic  effects  of 
ferroelectric  materials  having  the  perovskite  crystal  structure.  The 
model  describes  the  material  in  terms  of  its  crystal  lattice  dimensions, 
relative  positions  of  its  ions  within  the  lattice,  internal  electric 
fields,  and  ion  polarizabilities.  The  model  performs  a  self-consistent 
computer  solution  of  the  Lorentz  formula  for  the  crystal's  electric 
fields  and  polarizabilities  and  uses  the  results  in  solving  for  those 
components  of  the  dielectric  tensor  that  lie  along  the  crystal's  axes. 

Models  of  this  type  have  been  used  by  other  investigators  in 
attempts  to  calculate  refractive  indices,  dielectric  constants,  and 
other  basic  electromagnetic  material  parameters  (Ref  9»  14,  25).  These 
efforts  along  with  those  using  other  models,  such  as  the  thermodynamic 
treatment  of  Devonshire  (Ref  19:68)  have  met  with  mixed  success,  failing 
in  many  cases  to  calculate  correctly  basic  crystal  parameters  such  as 
the  natural  birefringence  in  barium  titanate  (BaTiO^)  and  lead  titanate 
(PbTiO^)  .  These  materials  are  optically  negative,  whereas,  calculations 
have  shown  them  to  be  optically  positive  (Ref  11:349). 

In  his  application  to  barium  titanate,  Wunsch  used  measured  lattice 
parameters  and  varied  slightly  the  approximate  polarizabilities  reported 
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in  the  literature  to  match  the  dielectric  constants,  refractive  indices, 
and  birefringence  in  a  free-field  environment.  Then  using  the  same 
basic  parameters  in  a  non-zero  field,  he  successfully  calculated  the 
linear  and  quadratic  electro-optic  constants. 

His  success  was  due  to  two  major  differences  incorporated  in  his 
model.  The 'first  and  most  important  was  to  use  an  anisotropic  field- 
dependent  form  of  the  polarizabilities.  This  form  was  obtained  by  a 
first  order  perturbation  quantum  mechanical  calculation  with  the  final 
result  being  expressed  in  terms  of  the  isotropic  zero  -field  polariza¬ 
bilities.  The  second  change  was  to  depart  from  a  pure  ionic  treatment 
by  incorporating  a  numerical  factor  that  simulates  the  degree  of  covlency 
of  bonds  between  ions.  This  factor  essentially  reduces  the  charge  at 
each  ion  site  by  the  amount  of  electronic  charge  that  is  shared  with  the 
oppositely  charged  ion  in  bonding. 

The  interest  in  being  able  to  theoretically  determine  the  electro¬ 
optic  behavior  of  linear  or  non-linear  optical  materials  is  due  to  their 
extensive  use  in  laser  applications.  The  most  common  of  these  uses 
include  laser  switching  and  amplitude  or  frequency  modulation,  with 
others  being  investigated,  such  as  information  storage  and  visual  display. 

Thesis  Problem 

The  problem  addressed  here  is  to  take  the  point-dipole  model  that 
Wunsch  applied  to  BaTiO^  and  test  it  further  by  applying  it  to  PbTiO^  . 
PbTiOj  appears  to  be  an  ideal  material  for  further  tests  on  the  model 
for  the  following  reasons:  l)  PbTiO^  and  BaTiOj  have  been  studied  more 
than  other  perovskites  because  only  these  two  materials,  of  the  many 
belonging  to  the  perovskite  ferroelectrics,  have  the  tetragonal  crystal 
structure  at  room  temperature.  Although  there  are  not  as  many  experi- 


2 


GEp/PH/72-19 


mental  measurements  reported  on  PbTiO^  as  on  BaTiO^  ,  due  to  greater 
difficulty  in  growing  crystals  of  it  and  in  making  measurements  on  it, 
there  is  sufficient  data  on  it  to  make  many  of  the  necessary  comparisons 
with  the  model’s  calculations.  2)  The  crystal  structure  characteristics 
that  are  chiefly  responsible  for  the  ferroelectricity  exhibited  are  more 
pronounced  in  PbliO^  .  Any  failings  or  weaknesses  of  the  model  should 
therefore  be  more  easily  detected.  And,  3)  Although  at  room  temperature 
the  spontaneous  polarization  of  PbTiO^  is  approximately  four  times  larger 
than  that  of  BaxiO,  ,  its  natural  birefringence  is  approximately  three 
times  smaller  (Ref  19:95)  .  This  is  contrary  to  the  theories  relating 
these  two  properties.  ‘The  behavior  of  the  natural  birefringence  with 
temperature  is  also  quite  different  for  PbTiO^  .  Any  success  by  the 
model  in  predicting  these  as  yet  unexplained  behavior  (Ref  19:96)  can 
be  considered  as  strong  support  for  the  model. 

The  goals  to  be  accomplished  and  the  approach  to  be  used  may  be 
summarized  in  the  following  way:  Using  experimental  crystal  structure 
parameters,  approximate  ionic  parameters  are  varied  to  match  as  closely 
as  possible  the  experimental  refractive  indices,  dielectric  constants, 
and  spontaneous  polarizations.  Then  keeping  the  crystal  structure  and 
ionic  parameters  fixed,  an  external  electric  field  is  introduced  into 
the  calculations  and  the  new  values  of  refractive  indices,  dielectric 
constants,  and  polarizations  are  obtained.  The  additional  anisotropy 
of  these  values  over  the  zero-field  values  are  then  used  to  determine 
the  appropriate  electro-optic  constants. 

An  additional  goal  was  preparation  of  existing  computer  codes, 
capable  of  performing  the  above  calculations,  to  ran  on  the  Air  Force 
Weapons  Laboratory's  computer  system. 
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II.  Discussion  of  Model 


Basic  Equations 

As  mentioned  in  the  introduction,  the  model  that  Wunsch  extended 
in  his  dissertation  to  do  electro-optic  calculations  is  a  point-dipole 
model.  The  point-dipole  model  requires  that  the  crystal  lattice  dimen¬ 
sions  and  positions  of  the  ions  within  that  lattice  be  well  known.  It 
assumes  that:  1 )  The  ions  have  constant  electric  moments  equal  to  the 
statistical  expectation  value,  2)  The  field-polarized  ions  interact  as 
point-dipoles,  and  3)  An  individual  ionic  and  electronic  polarizability 
can  be  assigned  to  each  type  of  ion.  The  basic  problem  faced  by  the 
model  is  a  solution  for  dielectric  constants  along  the  crystal's  prin¬ 
cipal  axes. 

The  i-th  component  of  the  dielectric  constant  for  an  anisotropic 
material  in  given  by 

f„Ei.=  Ei+  «jr7pf  (1) 

b 

_  b 

where  K.  is  the  external  electric  field  in  the  i-th  direction  and  P. 

I  I 

is  the  polarization  in  the  i-th  direction  of  each  type  of  ion.  The 
polarization  of  an  ion  type  per  unit  volume  is  given  by 


.  ,b  b  „b 

N  a  F: 


(2) 


where  [\J  is  the  number  of  ions  of  each  type  per  unit  volume.  Similarly, 

_b  b 

|-  .  is  the  local  field  at  each  ion  and  or  is  the  polarizability  of  the 


1 

ion  in  the  i-th  direction. 
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The  local  field  at  each  ion  is  given  by  the  Lorentz  formula.  For 
a  non-cubic  material  composed  of  different  ions,  this  formula  is 


e.  + 

'  bj  *J  J  J 


where:  p*3  =  local  field  at  the  a-th  ion  in  the  i-th  direction 

Ej  “  applied  fiel^  in  the  i-th  direction 
|\  =  number  of  type  *  b*  ions  per  unit  volume 

,  ab 

I. .  =  lorentz  factors 

y 

OT  =  field-dependent  polarizability  of  a  type  *b'  ion 

J 


Substituting  for  the  polarization  in  equation  (l)  from  equation  (2) 
and  solving  for  the  dielectric  constant  along  the  principal  axes,  gives 


b  r, 


where,  for  only  one  ion  of  each  type  per  unit  cell,  f\j:=  JL.  where  A  is 

pb  A 

the  volume  per  unit  cell.  The  -----  (local  field  factors)  are  found  by 

Ej 

solving  equation  (3)  .  Since  the  polarizabilities  are  defined  as  being 
field-dependent,  a  closed-form  solution  is  difficult  to  obtain.  A  major 
goal  of  the  model,  then,  is  "  '--’lf-consistent  solution  of  this  equaiion 
by  computer  techniques. 
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For  a  cubic  lattice,  each  of  the  lorentz  factors  is  equal  to  i/3  . 
Therefore,  fcr  a  cubic  lattice  Kith  one  atom  of  each  type  per  unit  cell 
equation  (3)  can  be  written  as 


Assuming  very  small  local  fields  so  that  the  polarizabilities  can  be 
taken  as  constants,  equation  (3)  can  be  solved  for  the  local  field 
factors  as 


b 


Than  substituting  for  the  local  factors  in  equation  (4),  the  dielectric 
constant  can  be  written  as 


For  a  material  that  docs  not  have  a  simple  cubic  lattice,  the 
solution  for  any  component  of  the  dielectric  constant  differs  from 
equation  (7)  in  that  the  denominator  of  the  local  field  factor  does 
not  involve  a  simple  prodi-'t  of  1/3  times  the  sum  of  polarizabilities. 
It  involves  a  sum  of  the  products  of  the.  proper  lorentz  factors  and 


polarizabilities. 
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The  methods  for  calculating  the  lorentz  factors,  field-dependent 
polarizabilities,  coulomb  fields  and  other  required  parameters  are 
treated  in  the  following  sections. 

Polarizabilities 

The  polarizabilities  of  equation  (3)  are  defined  to  be  field- 
dependent,  however,  only  a  portion  of  the  polarizability  is  given  a 
field  dependence.  The  polarizability  is  separated  into  an  electronic 
and  an  ionic  part  as  per  basic  assumption  3)  of  the  point-dipole  model. 
The  electronic  polarizability  is  defined  as  the  tendency  of  the  elec¬ 
tronic  charge  cloud  arpund  an  ion  to  be  distorted  by  an  electric  field. 
The  ionic  polarizability  is  the  tendency  of  an  ion  in  a  crystal  to  be 
displaced  from  its  normal  position  with  respect  to  other  ions  of 
different  charge  resulting  in  the  formation  of  dipoles  (Ref  12:164)  . 

The  ionic  polarizability  is  not  an  experimentally  measureable  value. 

It  is  determined  instead  by  taking  the  difference  between  the  total  and 
the  electronic  polarizability.  The  total  polarizability  is  defined  and 
determined  by  the  Clausius-Ilossotti  equation.  This  equation,  which  may 
be  obtained  by  rewriting  equation  (7),  is 


In  this  equation  is  the  total  polarizability  (response  to  a  static 
or  a  low  frequency  alternating  electric  field).  For  alternating  electric 
fields  of  optical  frequencies,  the  response  would  come  from  the  elec¬ 
tronic  charge  distributions  of  the  ions,  and  at  these  frequencies, 
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equation  (8)  nay  be  written  as  the  Lorenz-Lorentz  equation 


where  the  square  of  the  refractive  index  is  also  known  as  the  optical 
dielectric  constant,  and  the  oce  are  the  electronic  polarizabilities. 

The  electronic  polarizabilities  that  are  obtained  through  equation 
(9)  are  isotropic.  The  anisotropic  field-dependent  polarizabilities 
required  by  the  model  are  arrived  at  through  a  quantum-mechanical  pertur¬ 
bation  calculation.  The  resultant  relation  is 


where 


(10) 


(11) 


fi  = 


S/2  ^2  2 

a0  m  F 


In  using  equation  (10),  OCq  is  replaced  by  the  isotropic  electronic 
polarizability  determined  by  the  Lorenz-Lorentz  equation. 
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Lorentz  Factors 

The  use  of  correct  lorentz  factors  for  non-cubic  crystals  is  a  basic 
distinguishing  characteristic  of  the  rrodel,  and  although  a  rigorous 
derivation  is  very  long  (Hef  27:155)  ,  some  discussion  is  required  to  gain 
an  appreciation  for  what  the  model  involves.  This  discussion  tries  only 
to  show,  as  briefly  as  possible,  the  conditions  under  which  the  fora  of 
the  lorentz  factors  is  determined. 

Consider  a  neutral  crystal  with  a  periodic  charge  distribution  that 
can  be  expressed  in  a  fourier  series 

P  =  JpUG,q  r  (.3) 

L 

and  in  terms  of  the  three-dimensional  Dirac  delta  function  as 

P  =I>k<5(r-rk] 

k 

The  electric  potential  of  this  crystal  may  also  be  written  in  a 
fourier  series 

<()  =  CL  e'qL’?  (is) 

L 


or  given  by  Poisson’s  equation  as 

V?0  =  -47 rP 


06) 
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Having  solver  for  the  fouxier  coefficient  p  of  equation  (13)  by 
using  the  expression  for  p  from  equation  (14),  the  fourier  coefficient 
Q  ^  of  equation  (15)  nay  he  determined  by  substituting  the  expressions 
for  P  ,  equation  (13),  and  for  (j)  ,  equation  (15),  into  equation  (16). 
Having  done  so,  equation  (15)  may  be  written  as 


or 


where 


l?rVVn 
A  LL  k  (qL  f 

L  k  1 


0  =  lek0(*-*kl 

k 


£7rV  e'q 

A  A-  (sL)2 


(17) 


(IS) 


(19) 


A  system  of  dipoles  may  be  formed  by  placing  a  charge  'e*  at  the 
origin  and  an  equal  opposite  charge  very  near  it  at  p  and  repeating 
the  arrangement  within  each  unit  cell.  The  potential  for  this  lattice 
of  dipoles  is  then  given  by  equation  (18)  as 


0  =  e^(r]~  e^{f-<5?  j 

.«  e^(r)-e^(r)+ 

*  e<5  n  M\f/  =  -M  V^  (20) 


where  M=-®(5f  is  the  dipole  moment. 
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Cotfteb  Fields 

The  electric  fields  (coulcssb)  existing  within  the  crystal  at  the 
different  ion  sites,  due  to  the  unsyssetrical  distribution  of  the  ions, 
are  calculated  in  e  sanser  sisular  to  the  calculation  for  the  crystal, 
dipole  fields  discussed  in  the  preceding  section.  Once  the  electric 
potential  0  for  the  lattice  of  charges  is  defined  by  equation  (18), 
the  field  is  given  by 

E  = -V0  (Z8) 


Pbr  a  non-cubic  non-syrssetric  lattice,  the  lattice  is  broken  down 
into  a  number  of  sue- lattices,  each  of  which  is  cubic  and  seets  the 
requirements  of  the  model's  development.  For  such  a  lattice,  the  field 
at  an  ion  site  is  found  by  adding  the  field  contributions  from  each  of 
the  different  sub-lattices. 


Spontaneous  Polarization 

The  spontaneous  polarization  is  equal  to  the  sun  of  the  electronic 
and  the  ionic  polarizations.  'The  electronic  contribution  by  each  ion  is 
due  to  the  deformation  of  its  electronic  charge  distribution  by  the  local 


field  and  is  given  by 


pa_  aaFa 

p  ~  ~K~ 


where  or  is  the  anisotropic  f 5. eld-d ep end ent  electronic  polarizability. 
The  ionic  polarizability  has  been  defined  within  a  crystal  (page  7) 
in  terns  of  the  displacement  of  an  ion  with  respect  to  other  ions  in  the 
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lattice  The  ionic  sposstanesicss  polarization  is  therefore  calculated 
in  terns  of  the  ion  shifts  azay  from  their  symmetrical  lattice  poult  lens. 
The  equation  for  calculating  the  polarization  of  each  ion  type  per  unit 
cell,  in  tee  direction  of  its  displacement  Z  ,  is 


p  -  Nqzd 
*”  A 


where:  ^  *  number  of  ions  per  unit  cell 

*  the  electronic  charge 
2  *  degree  of  ionization  of  the  ion 
d  -  the  ion's  shift  froa  its  symmetrical  cubic  position 

Electro-optic  Constants 

When  the  perovskite  crystal  is  subjected  to  an  external  electric 
field,  the  electronic  charge  distribution  of  the  individual  ions  suffer 
a  deformation  in  addition  to  any  that  nay  already  exist  due  to  internal 
coulomb  fields.  This  defoliation  results  in  a  change  in  the  electronic 
polarizability  in  the  direction  of  the  applied  field.  In  a  sinner 
identical  xith  the  natural  birefringence,  the  induced  anisotropy  in 
electronic  polarizabilities  results  in  an  induced  birefringence.  This 
Induced  birefringence  can  be  expressed  as  a  sun  of  a  linear  and  quadratic 
tera  in  the  electric  field  strength,  with  the  terms  being  called  the 
Pockel  and  the  Kerr  effect,  respectively. 

In  the  crystal's  cubic  phase,  the  linear  tern  is  absent  due  to  the 
crystal's  symmetry  (Ref  29:301).  The  constant  for  the  quadratic  effect 
can  be  expressed  (Appendix  A)  as 

d  o  —  2 An 

"ll  ^12  “  n3p2  (31) 

O  LZ 
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where  An  is  the  induced  birefringence  and  no  is  the  zero-field  refrac¬ 
tive  index  along  either  of  the  principal  axes.  The  above  constant  is  also 
expressed  in  terns  of  the  induced  polarization,  rather  than  the  applied 

f 

electric  field.  This  relation  is  given  by 


q  _  q  s  _2.An 

yil  y!2  n3p2 

0*2 


The  induced  polarization  that  goes  into  the  above  relation  is  calcu¬ 
lated  by  equation  (29)*  where  the  polarizability  is  now  a  sum  of  both 
the  electronic  and  the  ionic  polarizability,  and  the  field  is  the  local 
field  solution  of  equation  (3)  when  an  external  electric  field  is  present. 

In  the  tetragonal  phase,  the  linear  effect  is  dominant  and  now  the 
quadratic  contribution  is  neglected  to  write  the  expression  for  the 
linear  electro-optic  constant  (Appendix  A)  as 


r  is 


^3  = 


z[ ne-  "z  1 


"eEz 


In  this  expression,  no  and  ne  are  the  zero-field  refractive  Indices 
along  the  'a*  axis  and  *c*  axis,  respectively,  is  the  refractive  index 
along  the  *c*  axis  in  the  presence  of  the  applied  field. 
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HI.  Perovskite  Ferroelectrics 


General  Properties 

f  The  perovskite  ferroelectrics  are  a  large  number  of  materials  having 
the  chemical  formula  ABO^  ,  where  generally  A  is  a  metal,  B  is  a  transition 
element,  and  0  is  oxygen.  This  family  is  a  member  of  still  a  larger  family 
entitled  the  oxygen-octahedra  family,  and  is  itself  divided  into  smr'ler 
families.  Barium  titanate  and  lead  titanate  belong  to  one  of  the  smaller 


families  where  the  A  atom  is  bi- valent  and  the  B  atom  is  tetra- valent. 


The  perovskite  ferroelectrics  meet  the  additional  requirement  that  their 
structures  are  similar  to  that  of  perovskite  (CaTiO^)  at  room  temperature 
(Figure  l).  This  material  from  which  the  family  took  its  name  has  since 
been  shown  to  have  a  slightly  distorted  cubic  structure  at  room  temperature 
(Ref  15:10). 


Fig.  1 .  Cubic  Perovskite  Structure 


The  materials  must  of  course  also  satisfy  the  definition  of  a  ferro¬ 
electric.  Although  at  first  this  might  seem  a  trivial  matter,  it  is  not, 
as  has  been  shown  by  a  number  of  attempts  to  clearly  define  the  properties 
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of  ferroelectrics.  This  problem  continues  even  up  to  the  present  due 
to  the  difficulty  of  classifying  the  many  new  solutions  and  mixtures  of 
these  materials.  However,  the  following  definition  discussed  in  length 

t 

by  Bogdanov  (Ref  2*591 )  appears  to  be  the  most  basic  and  clearest  - 
"ferroelectric  materials  are  those  which  in  some  range  of  temperature 
show  spontaneous  polarization,  the  direction  of  which  may  be  changed  by 
an  external  electric  field." 

The  perovskites  do  not  generally  have  the  cubic  structure  of  Figure  1, 
but  their  structures  can  be  described  in  terms  of  the  cubic  structure  as 
slight  expansions  or  contractions  along  an  axis  or  diagonal.  For  barium 
titanate,  these  structures  are  tetragonal,  orthorhombic,  and  trigonal 
over  different  temperature  ranges.  The  tetragonal  structure  is  the 
simplest  and  consists. of  an  expansion  along  one  axis  ('c'  axis)  and  a 
contraction  along  the  other  two  axes  ('a'  and  'b'  axes).  At  high  enough 
temperatures,  these  ferroelectrics  do  attain  a  cubic  structure  and  lose 
their  ferroelectric  properties.  This  transition  temperature  is  called 
the  'curie'  temperature  and  the  higher  temperature  phase  is  called  the 
paraelectric  phase  in  analogy  with  the  ferromagnetic  to  paramagnetic 
transition. 

In  applying  the  point-dipole  model  to  the  perovskite  crystal,  the 
crystal  lattice  must  be  divided  into  sub- lattices  which  have  the  perio¬ 
dicity  assumed  in  the  development  of  the  model.  The  sub-lattices 
required  for  the  perovskite  crystal  are  (see  Figure  l)t 

1.  Barium  ion  lattice 

2.  Titanium  ion  lattice 

3*  Oxygen  I  ion  lattice  (in  yz  plane) 

4,  Oxygen  II  ion  lattice  (in  xz  plane) 

5.  Oxygen  III  ion  lattice  (in  xy  plane) 
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Lead  fi Senate 

Lead  t&taaste  exists  in  the  cubic  perovskite  structure  above  450*0 
ard  in  a  tetragonal  structure  below  this  temperature.  This  high  curie 
teaperature  has  been  attributed  to  a  higher  polarizability  of  the  lead 
ion  than  that  of  the  baritn  ion  (Ref  1:1433).  Although  sone  anaaolies 
were  seen  in  dielectric  constant  measurements  at  tenperatures  lower 
than  roan  temperature  (Eef  13)»  it  has  not  been  established  that  it 
goes  through  any  additional  transitions  at  lover  tenperatures. 

Hie  transition  to  the  tetragonal  structure  involves  acre  than  the 
simple  expansion  and  contraction  along  +ue  crystallographic  axes  of  the 
cubic  structure  as  mentioned  above,  it  also  involves  shifts  cf  the  ions 
away  from  their  symmetrical  lattice  positions.  Bo  satisfactory  explan¬ 
ation  exists  for  the  condition  causing  these  shifts,  but  the  situation 
does  result  in  internal  fields  that  are  no  longer  zero  at  the  various 
ion  sites.  The  fields  at  the  ion  sites  tend  to  polarize  the  ions,  pro¬ 
ducing  dipoles  whose  fields  further  affect  the  fields  and  ions  at  other 
lattice  sites.  Since  the  polarizabilities  depend  on  the  crystal's 
electric  fields,  it  is  these  deviations  from  a  cubic  perovskite  lattice 
that  can  be  related  to  the  ferroelectric  properties  of  the  barium 
titanate-type  perovskites. 

The  basic  PbTiO^  crystal  dimensions  were  determined  by  Shirane  and 
Hoshino  (Ref  22:266)  in  1950  by  x-ray  diffraction  analysis.  Their  results 
for  crystal  lattice  dimensions  versus  temperature  in  the  tetragonal  and 
cubic  phase  are  shown  in  Figure  2  ,  and  are  listed  along  with  other 
crystal  physical  properties  in  Table  I  . 

The  ion  positions  have  been  determined  by  Shirane  and  co-workers 
(ref  21:136)  and  the  accepted  model  is  shown  in  Figure  3  . 
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Table  I. 


PbTiO^  Crystal  Lattice  Parameters 


Temperature 

(°c) 

c  axis 

a-) 

a  axis 
(A) 

c/a 

Unit  Cell  Vol. 
(?) 

30 

4.140 

3.694 

1.063 

62.77 

110 

4.125 

3.899  .  • 

1 .058  .  • 

62.72 

200 

4.104 

.  3.907 

1.050 

62.65 

300 

4.081 

•  3.918 

1 .041 

62.65 

325 

4.073 

3.921 

1 .039 

62.61 

350 

4.066 

3.923 

1.03V 

62.58 

400 

4.046 

3.930 

.  1.029 

62.49 

410 

4.041 . 

3.933 

1.028 

62.50 

425 

4.036 

3.934  ’ 

1.026 

62.45 

440 

4.030 

3.936 

1.024 

62.43  > 

450 

4.028 

3.937  , 

1.023 

62.44 

460 

4.021 

3.938 

1 .021 

:  62.37 

470 

4.016 

'  3.942 

•  1,019 

62.31 

480 

■  4.011 

3.942 

1.018 

62.31 

490 

3.958 

3.958. 

1.000 

62.03 

510 

3.959  , 

3.959 

1.000 

62.05' 

525 

•  3.960 

'>.960 

1.000  V 

62.09 

535 

3.960 

5,960 

'  1.000 

62.10 

*  800  ■ 

3.964 

964 

1.000 

1 

*  (Ref  4:65) 

_ L. 

(Ref  22; 26?) 
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<5  zTi  “  +  0*040 


<5v "  <57o  "  6z o  “ 

i  i:  hi 

(Note*  (5z*s  are  ratio  of  the  atom’s 
shift  distance  to  ’c’  axis  dimension) 


Fig.  3*  PbTiO^  Tetragonal  Crystal  Lattice 
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A  very  distinguishing  characteristic  of  PbTiO^  among  the  perovsldte 
ferroelectrics  is  the  natural  birefringence  exhibited  by  the  tetragonal 
phase.  This  natural  birefringence  has  been  measured  fcv  Shirane  (Eef  21 : 

i 

140)  and  by  Fesenko  (Ref  5:55).  Figure  4  shows  the  results  obtained  by 
Shirane  for  FbTiO^  .  To  emphasize  the  distinguishing  behavior.  Figure 
5  shows  the  natural  birefringence  exhibited  by  BaTiO^  (Ref  19:95), 
which  is  more  representive  of  other  perovskite  ferroelectrics.  The 
sharp  difference  in  values  at  room  temperature  and  in  the  behavior  with 
temperature  can  be  easily  seen. 

It  must  be  pointed  out  that  it  is  more  difficult  to  work  experi¬ 
mentally  with  FbTiO^  than  with  BaTiO^  ,  due  to  greater  difficulty  in 

producing  single-domain,  or  even  multi-domain,  single  crystals  of  FoTiO_  -. 

5 

Measurements  are  hampered  by  the  higher  fields  required  to  reverse  the 
polarization  in  PbTiO^  and  by  conductivities  large  enough  to  influence 
dielectric  hysterises  measurements.  Because  of  these  difficulties, 
some  of  the  reported  parameter  values,  needed  for  comparisons  with  the 
model’s  predictions,  have  involved  approximations  that  reduce  their 
credibility.  This  is  true  for  the  electro-optic  constants.  Other 
parameters,  such  as  the  natural  birefringence  and  other  basic  properties 
are  well  established,  however,  and  should  provide  a  good  test  for  the 


model. 
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Fig.  4.  Natural  Birefringence  in  Tetragonal  PbTiO_ 
(Experimental) 


Fig.  5. 


Natural  Birefringence  of  Tetragonal  BaTiO 
(Experimental) 
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IV.  Lorentz  Factors 

As  mentioned  in  discussing  th®  model,  the  lorentz  factors  for  the 
perovskite  crystal  cannot  he  approximated  by  a  value  of  i/3-  The 
lorentz  factors  mast  he  calculated  with  every  lattice  change,  with  the 
exception  of  changes  in  the  cubic  phase.  Moire  specifically,  in  the 
tetragonal  phase,  where  every  change  in  temperature  produces  a  change 
in  the  c/a  ratio,  new  lorentz  factors  mist  be  calculated  at  each  new 
temperature.  In  the  cubic  phase,  the  three  crystal  axes  dimensions  are 
equal  and  increase  equally  with  temperature,  therefore,  once  the  lorentz 
factors  have  been  calculated  at  any  temperature  within  the  cubic  phase, 
these  sane  factors  can  be  used  at  any  temperature  in  that  phase. 

The  crystal  dimensions  and  ion  shifts,  used  in  calculating  lorentz 
factors  at  the  different  temperatures  considered,  were  taken  from 
Ihble  I  and  Figure  3  ,  respectively.  The  one  exception  was  the  value 
of  the  *a*  axis  at  room  temperature,  where  a  value  of  3-899  A  was  used 
as  an  average  of  the  values  measured  by  various  investigators. 

The  lorentz  factors  for  the  tetragonal  phase  at  rooa  temperature 
and  for  the  cubic  phase  at  the  curie  temperature  are  given  in  Table  II 
and  Table  III,  respectively.  The  cross-terms  are  not  Included  because 
they  are  equal  to  zero  due  to  the  crytal  lattice  symmetry. 


23 


GEP/PH/72-19 


i 


owa^w 

g*o((\«0 

©  t-\o  ct  n 

ffk 

OotrtO'C 
U>tC\e>  {V  © 
©  ©  WJfVJ-A 

•  •  •  •  • 

ooooo 


IT  IflOVO  © 

in^o  r  i 

1*  ©  OlfAp. 

^0»- 

UW^CtpJ 

cn  in  -c-m  o. 
o  -tr  •«#■  crvvo 
rvKvifvtf  o 

lAh^toni 

c-ot'-mm 

•  •  •  •  a 

o  o  o  o  o 

I  I 


Ift^tvo  jvo 

in  CM  K\  03  V- 
OOH\I*-VO 
■w-  CM  03  (r\ 

tn  ©  m  -<^  © 
tDOO'.'tlft 
tn  -Vd*  in  LTV 

e»  *-  ©  «-  co 

tn  f~-  tn  c~-  in 

•  •  •  •  • 

owooo 


cMvo-^m*- 

03  (Tl  -vp-  -<£-  G\ 
©  m  CM  ©  -r- 

O'ojcnow 

f  (J>(V  ■<*■«> 

•"-K'xcoirxc— 
3-  cn  o  in 
00  ^'tirvvo 
O  vo  ■«-  n-  in 
m  tn  t-  co  co 

•  •  •  1  • 

O  O  CM  o  o 

I  I 


H 


CO 

u 

o 

-p 

o 

ui 

N 

-P 

c 

d> 

h 

o 

PI 


©  cm  m 
cn  03  m 
tn  ©  in 
CM  03  O 

Oil-  T- 

tn  •»-  in 

03  t-  co 

M-OOm 
10  O  Is- 
tn  tn  tn 


03  in 
in  r— 
•»d-  vo 
«-  o 
in  co 

03  © 

o  -a- 
tn  © 
in  in 
t-  vo 


OOOOO 

l 


*H  M  M 
W  H 
H 

£>  -H  I?  U  X 
O  O  O 


£ 


a j 

u 

o 

-p 

o 

o 

&. 

N 

-P 

C 

d> 


in«-oo<s 

•»—  or»e-»a  cv 

W  ^  >  r-  ri 

O  w^-vow 

scsmsv 

©  ?-  Cu3  sn 
m-uvvo  «n  o 
cvootriv- 
«r»  m«M  ©  vo 
<0  ©  !n  insn 

•  •  •  •  m 

OOOOO 


snM-ovova 
in-c-vo  ewe 
in  cm  c-.yv^ 
O  03  CM  VO 
r-W  M-  C73  m 
in©  -vj- m vo 
©  O  v  C3  tn 
x\  -a-  m  -*Mn 

f**-T-UM3 
mt^t-min 
•  •  •  «  • 
owooo 


03  mvo  o  © 
in -a-  cm©  t>- 
Mtoncnh 
»-OCM*-M- 

in  03  •sj-  cm 

OV'A>3ttCft 

OMCftM-O 
IOK\<-!ftO 
in  r-  vo  t—  cm 
r-  ©  K\t-K3 

•  •  •  •  • 

OOOOO 

I  I 


CM  vo  in  «— 

03  03  M-  \1  03 
VO!<3VOW»- 
C3  CM  O  G3  CM 
03  -vf  CM  © 

T-niniDP- 
*-  C3  -si-  o  in 

©  -vl-  K3  -vl-  VO 

o  vo  r—  «—  in 
mm©  t-© 
•  •  •  •  • 
O  O  O  CM  o 

I  I 


vo  cm  03  in  in 
03  03intn»- 
tn  ©  -a-  m  © 
(M  03  1-  O  O 
03  ^  in  © 
tn  -r-  C3  in  © 
03  »-  o  ©  -a- 

a-  ©  tn  tn  © 
vo  o  in  c—  in 
tn  tn  i—  tn  © 


OOOOO 

l 


M  M  H 
M  M 
H 

.0  -H  x  i?  ^ 

&■  EH  O  O  O 


N 


0} 

f-S 

O 

-P 

o 

fi 

N 

■P 

a 

<d 

it 

o 

PI 


owm«n© 

CVS  €3  ■»-  •»-  O 
CM  nvoifl  CM 

^vCOEtUV 

vom  ©  ©  «- 

g^mCM  CM  CM 

vojcmtjvo 

jn  n- 1~-  cm 


owooo 


msioon 

O'  G3  CM  o  ■*- 
O  m  CA  CM  VO 

03  ©  m  in  o 
in©  ©,-  © 
m  cm  o  CM  CM 
I>-VO  CD 

ooootn 

CM  sdK\?'J<3 


VO 


1  mw  e» 


OOOOO 


in  03 

03  03 

o  in 

C3© 

in© 
in  cm 
r-vo 

o  o 
cm  -a- 
u>  © 


raotn 
o  w  t- 
CM  03  VO 

in  moo 

-r-  ©  © 
CM  ©  (M 
-^-  ©  ,- 
O  O  in 
tn  in 
cm  -a-  r~- 


OOOOO 


VO  ©  03  03  CM 
«-  O  03  03© 

vo  cm  in  in  m 
ointooM- 

©  «-  ©  ©  CM 
CM  CM  CM  in 
c—  »-  ©  ©  .- 
•n  poo  tn 

tn  cm  ©  ©  t'- 
•  •  •  •  • 

O  O  O  O  CM 

I  I 


©  © 
0  * — 

CM  © 

in  o 

T-  © 

CM  t> 
*—  C 

O  tn 
t~-  © 
cm  tn 


tn  in  tn 
03  03  CM 
o  o  CM 
03  03  «- 

in  in  © 
in  in  t- 

C-  C-  03 

00© 

N  Wt- 

©  ©  tn 


OOOOO 

1 


M  M  M 
M  M 
M 


X)  -H  £?  i? 

CM  EH  O  O  O 


24 


r~ 


CEP/Fg/72-«S 


H 

H 


0 

t-t 

"S 

e< 


s 

p 


"8 


o 

nrl 

o 

ft 

£ 


•2  **" 
tr 


h 

o 

•p 

o 

& 

fa 

N 

P 

C 

ffl 

P« 


H 

•>««-«  rv 

cm  e-*-  —  m 

1-1 

nir^-^rv 

nnMcvn 

*e-40  cm  cm  m 

H 

ao  o  *f  m 

m«Mv©  mm 

vo  -emmm 

« 

rl<V«r«\ 

cm  vo  *-  *-  m 

© 

r-  ©»*«!  rv 

*-  oca  *-  an 

cm  o  *-  *-  m 

22 

W  atM^rv 

NH4««m 

*r  *-  cm  cm  m 

crtotan 

cd  m  40  cs  m' 

40  C—  CD  CD  m 

p-wr-wn 

t-vo  me- an 

mcM  r- 1-  m 

4©  OOr.»\ 

vo  cd  m*©  m 

m  t-4©  40  m 

•  •  •  •  • 

•  •  •  •  • 

•  •  «  •  • 

ooooo 

©CJOOO 

©CMOOO 

*-  'M'*r-mCM 

cm  o»T-r\»- 

r-MCMm*- 

H 

cu  mcM  m-**- 

4t4C  CM  mCM 

cMmvtmcM 

H 

♦  O^K\© 

fflOviWM 

Moonvr 

•< 

n  cm  iorvw 

40  M-rtnrv 

m  cm  vo  mm 

C 

r-rv^rvw 

CM  40  -r-  m*p- 

*-  mcM  m*~ 

CD 

0 

*-  ot^nrti 

cm  oo  *-  m*- 

*-  o»cm  m*- 

m 

te 

CM  ITVCM  KV  < 

*0-r-  CM  K\CM 

cm  in  *4- m  cm 

<T» 

h 

CD  m  30  m  v.'> 

vo  r-  co  m  o 

to  mvo  m  o 

• 

X 

in  CM  MOH 

c—  vo  inmt— 

m 

o 

40  CO  40  m  m 

me- vo  m4o 

vo  cd  mmvo 

•  0  0  •  • 

•  •  •  •  • 

***** 

H 

OOOOCjl 

O  CM  O  O  O 

ooooo 

(UO|l\r-f 

T-M-nr-  CM 

*-  -M-mCM  *- 

M-UnNOI 

cu  mmeu  *4- 

cm  mm-M-cM 

M 

CD  OK\  ■42-41- 

onvjn 

*4-0  moo  *»- 

40  M-KM«M<\ 

m  cm  m  mvo 

m  cm  mvo  m 

c 

CMVOm*-*- 

«-  mm*-  cm 

*-  m  mcM  *- 

e 

CM  CO  m*-  *- 

*-  O'  m  «—  cm 

*-  cnmni  *- 

*3 

cm  in  m  cm  ■41- 

cm  in  m  *4-  cm 

40  f—  K\  CD  CD 

go  m  m  co  40 

co  mmvo  co 

8  ^ 

X 

in  <vi  m>-  c- 

t-vome-in 

t-vo  mint— 

1 

o 

m  e-mvo  vo 

40  cd  mvo  m 

vo  co  mmvo 

8  1 

•  *  *  •  • 

•  •  •  *  • 

***** 

owooo 

OOOOO 

ooooo 

i 

K\  m  CO  *4-  *«t 

mnM-OVM- 

KVKVtvfH 

m  in  vo  m  m 

mmm4o  m 

mmmmvo 

g 

m  m  o  o  o 

mmo  o  o 

mmo  o  o 

8  10 

m  m  *4-  cm  cm 

m  m  cm  *4-  cm 

mmcM  cm  *4- 

1  cr» 

-ri 

mmvo  in m 

mm  mvo  m 

mmmmvo 

B 

C 

mmco  04  0* 

mmovoo  cn 

m  m  cr»  cv  co 

1  •' 

3 

m  m *-  in  tn 

mm  in  *-  in 

mminin*- 

-p 

m  inf-in  tn 

mmmt-m 

m  m  m  m  c- 

i  m 

-rl 

m  mc-vo  vo 

mmvo  cm  vo 

m mvo  vo  cm 

\  ■q 

t-t 

mmiM  co  co 

mmco  c-cd 

mmcocot- 

1  K1 

•  •  •  •  • 

•  •  •  •  • 

•  •  •  •  • 

O  O  CM  O  O 

O  O  O  CM  o 

o  o  O  o  CM 

« 

a 

1  1 

1  1 

1  I 

O 

m  m  cm  t-  r 

mm*-  cm  *- 

mm*-  *-  cm 

tn  tn  *4-  cm  cm 

m  m  cm  vt  cm 

m  mcM  cm  *4- 

to 

KMfirovtvj- 

mm*a-  co  *4- 

mm*4-*i-co 

c 

tn  in  40  m  tn 

m  m  tn  4o  m 

m  m  m  mvo 

<D 

TJ 

m  tn  cm  *-  *- 

mm*-cM*- 

m  m  *-  *-  cm 

a 

Gj 

Kin«t-r 

m  m  *-  cm  *- 

m  m  *-  *-  cm 

*H 

© 

KV  tn  -d-  CM  CM 

m  m  cm  *4-  cm 

mmcMCM*4- 

o 

hI 

tn  tn  40  co  co 

m  m  co  vo  co 

m  m  co  oo  vo 

m  tn  in  c-  r- 

m  m  c-  in  t- 

mmc-r-in 

© 

« 

tn  mtnvo  vo 

m  m  ’o  tn  vo 

14 

m  mvo  vo  m 

O 

X 

•  •  «  •  « 

>» 

«  »  *  •  • 

N 

•  •  •  •  • 

•ri 

1 

OOOOO 

ooooo 

ooooo 

P 

-P 

1 

1 

1 

1 

1 

d 

o 

CJ 

CJ 

Pi 

u 

o 

o 

o 

p 

-p 

p 

o 

CJ 

CJ 

d 

d 

as 

to 

to 

fa 

fa 

fa 

£3 

c 

o 

o 

N 

HUM 

N 

HHH 

N 

H  H  M 

-i  i 

•rl 

M  H 

P 

M  M 

P 

M  M 

c 

M 

C 

H 

PS 

M 

1 

•- 

© 

>»>!>> 

© 

© 

d 

Pi 

Pi 

P> 

Pi 

^J  -rl  X  ^ 

o 

r©  -H  X  X  X 

o 

.©  -rl  X  X  X 

o 

I 

h? 

fa  C-I  o  o  o 

1-3 

fa  EH  a  o  o 

h! 

fa  E-t  O  O  O 

25 


CEP/SS/72-19 


Y»  Sdtotoe  5ier  gad  Saectregic  Polarizabilities 

Electronic  polari  zsbi  li  ties  reported  is  the  literature  hare  bees 
determined  by  the  requirement  that  they  predict  correct  refractive  indices 
through  the  Loreas>Lsxesiz  equation.  She  procedure  is  to  consider  a 
mater  of  different  materials  containing  different  combinations  of  the 
ions  whose  polarizabilities  are  desired.  By  *  trial-and— error 1  ,  elec¬ 
tronic  polarizabilities  are  assigned  to  the  individual  ions  and  then 
sunned  tr,  satisfy  the  molecular  compositions  of  the  different  caterials. 

She  correct  polarizabilities  are  determined  to  be  those  that  are  cost 
consistent  in  predicting  the  correct  refractive  indices  for  the  naterials 
considered.  A  study  of  this  type  has  been  done  by  Tessnan,  Kahn,  and 
Schockley  (Ref  26)  for  a  nuaber  of  atons  including  lead,  titanium,  and 
oxygen.  Polarizabilities  obtained  by  this  method,  however,  can  be 
considered  to  be  only  approximately  correct  for  two  reasons.  First, 
the  Lorenz-Lorenr  ;  equation  has  been  shown  to  hold  for  materials  with 
a  cubic  lattice,  but  in  the  above  study,  racy  non-cubic  materials  were 
included.  And  second,  an  ion's  electronic  polarizability  should  not  be 
expected  to  be  a  constant,  it  should  depend  on  its  environment  (Ref  12:166), 
Slight  variations  of  the  reported  values  may  therefore  be  necessary  for 
good  predictions  of  the  refractive-  index  for  a  particular  material. 

The  approximate  polarizabilities  determined  by  the  above  method 
are  here  used  and  varied  in  equation  4  >  using  correct  lorent.z  factors, 
to  match  the  experimental  cubic  phase  refractive  index  at  the  curie 
temperature.  The  electronic  polarizabilities  of  oxygen  and  titanium 
that  satisfy  the  tc ' ractive  index  of  BaTiCL  have  been  determined  by 
Wunsch.  These  valuer,  are  2.322  P  and  0.19  P  ,  respectively  (Ref  27:30). 
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S&ac®  the  environments  of  there  loss  ore  iiesticsl  in  the  cubic  pease 
of  both  materials,  with  the  exception  of  replacing  barium  'ey  lead,  their 
electronic  gol  ari  zgfeili  ties  were  held  constant  at  the  above  values  and 
the  reported  lead  value  vas  varied  to  satisfy  a  refractive  index  equal 
to  2.72  (Jtef  24:4792)  for  cubic  TVSiO^  at  its  curie  temperature. 

Vesssan  (lef  26:893)  reports  the  electronic  polarizability  of  lead  to 
be  4.79  A3  ,  while  Fedulov  (Ref  4:63^  reports  a  value  of  4.32  X3  . 

In  doing  the  calculations,  the  lead  value  was  varied  fro*  below 

•3  «3 

3.00  A  to  5.00  A  •  The  first  set  of  calculations  done  for  a  lattice 

dimension  of  3.960  A3  shoved  that  the  lead  electronic  polarizability  had 

•3 

to  be  extended  below  3-0  A  before  the  2.72  refractive  index  value  could 
be  obtained  (Figure  6).  The  other  curves  in  the  figure  correspond  to 
larger  values  of  the  lattice  dimensions  and  give  an  idea  of  the  error 
that  nay  be  expected  in  the  calculated  refractive  index  due  to  errors 
in  the  measured  lattice  dimensions.  Since  the  measurement  error  is  less 
than  the  difference  between  the  lattice  dimension  corresponding  to  the 
upper  and  lower  curve  in  the  figure,  the  calculation  error  due  to  uncer¬ 
tainty  in  the  lattice  dimension  should  be  small.  The  lead  electronic 

polarizability  required  to  match  the  refractive  index  at  a  lattice 

•  »3 

dimension  of  3.958  A  is  2.90  A  .  Table  IV  summarizes  these  values,  which 
will  he  used  throughout  the  remaining  calculations,  both  in  the  cubic  and 
tetragonal  phase.  In  the  tetragonal  phase,  these  values  are  the  isotropic 

electronxcpolarizabilities  that  go  into  equation  10  . 

°3 

The  2.90  A  lead  value  is  much  lower  than  the  two  values  reported 

in  the  literature.  A  close  study  of  Tessman's  value,  however,  shows 

°3 

that  2.90  A  may  not  be  an  unreasonable  value.  The  original  source  of 

®3 

the  4.32  A  value  was  not  available  for  further  consideration. 


Fig.  6.  Refractive  Index  vs.  Lead  Electronic  Polarizability  in 
Cubic  PbTiOj  (Theoretical) 


28 


GEP/PH/72-19 


Table  IV. 


Electronic  Polarizabilities  for  Ions  in  PbTiO- 

5 


Ion 

Electronic  Polarizability 

Pb 

2.90  A3 

Ti 

0.19  I3 

0 

2.322  A3 

GEP/PH/72-1 9 


VI.  Dielectric  Constants  and  Total  Polarizabilities 

The  procedure  used  in  determining  the  total  polarizabilities  that 
have  been  reported  in  the  literature,  is  similar  to  the  one  for  electronic 

i 

polarizabilities  from  the  previous  section.  By  using  the  Clausius- 
Mossotti  equation,  total  polarizabilities  are  established  for  a  number 
of  ions  by  requiring  that  they  be  consistent  with  the  dielectric  constants 
of  a  number  of  materials  involving  the  ions.  This  has  been  done  by 
Roberts  (Ref  17,  18)  . 

Here,  as  with  the  electronic  polarizabilities,  total  polarizabilities 
are  established  by  equation  4  using  correct  lorentz  factors  to  match  the 
dielectric  constant  of  cubic  PbTiO^  at  the  curie  temperature  and  of 
tetragonal  PbTiO^  at  room  temperature.  The  experimental  dielectric 
constants  to  be  matched  are  10,000  for  the  cubic  phase  (Ref  16:41 )  and 
210  and  120  for  the  tetragonal  'a*  and  'c'  axis,  respectively  (Ref  7:549). 

As  with  the  electronic  polarizabilities,  Wunsch  (Ref  27:87)  has 
shown  that  total  polarizabilities  of  0.875  A  and  2.8518  A  for  titanium 
and  oxygen,  respectively,  satisfy  the  dielectric  constant  of  cubic 
BaTiO^  .  An  attempt,  however,  to  match  the  dielectric  constant  of  cubic 
PbTiO^  at  the  curie  temperature  by  using  the  above  values  for  titanium 
and  oxygen  and  varying  the  lead's  value  around  Roberts'  value,  7.1  A 
(Ref  17:866)  resulted  in  negative  dielectric  constants. 

The  negative  dielectric  constants  may  be  explained  by  considering 
equation  4  ,  where  the  presence  of  negative  local  field  factors  (summa¬ 
tion  term  in  equation  6,  being  larger  than  one)  can  result  in  negative 
dielectric  constants.  At  any  one  temperature,  the  lorentz  factors  and 
unit  cell  volume  are  fixed  and  therefore  only  the  polarizabilities  may 


30 


Qs&fssfiz-n 


be  varied  to  correct  the  situation.  The  loreniz  factors,  however, 

certain  both  positive  asd  negative  terns  and  so  it  is  not  clear  which 

polarizabilities  should  he  increased  or  decreased  to  produce  the  desired 
*• 

results.  It  is  necessary,  therefore,  to  solve  the  complete  set  of 
equations  for  every  variation,  and  this  requires  a  very  large  number 
of  calculations.  To  sake  this  possible,  the  various  computer  programs 
responsible  for  the  different  parts  of  the  model  were  combined  into  one 
program  with  the  necessary  modifications  and  additions  to  allow  the 
variations  to  be  done  efficiently. 

Once  this  was  done,  the  polarizability  variations  consisted  of  two 

\  ®3  °3 

parts:  1 )  vrrying  the  values  for  lead  from  4.5  A  to  7.7  A  and  for 

•3  *3 

oxygen  from  2.672  A  to  2.852  A  while  holding  the  titanium  value 
°3  \ 

constant  at  0.19  A  ,  and  2)  varying  the  lead  and  and  oxygen  over  the 
same  range  as  in  l),  but  for  values  of  titanium  polarizabilities  that 
were  increased  and  decreased  by  5,  10,  and  50^  of  the  ionic  polarizability. 

The  results  of  the  polarization  variations  were  not  as  expected  or 
hoped  for.  No  variation  was  found  that  could  simultaneously  predict 
the  experimental  cubic  and  tetragonal  dielectric  constants.  Instead, 
se^-ral  variations  were  found  to  sat.'sfy  the  cubic  dielectric  constant, 
and  although  only  one  region  of  polarizabilities  satisfied  simultaneously 
the  tetragonal  'a'  and  'c*  axes,  several  satisfied  the  'a'  axis  alone. 
Polarizabilities  in  the  region  thr.t  satisfied  both  tetragonal  axes, 
however,  resulted  in  nixed  positive  and  negative  local  field  factors, 
and  this  has  been  suggested  as  indicating  antiferroelectric  tendencies 
(Ref  28). 

The  promising  polarizability  combinations  were  used  in  extending 
the  dielectric  constant  calculations  with  respect  to  temperature  in 
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both  the  cubic  and  tetragonal  phase*  In  the  cubic  phase,  these  polari¬ 
zabilities  were  also  used  in  calculating  the  electro-optic  effects.  It 
was  hoped  that  these  calculations,  which  are  discussed  in  the  following 
sections,  would1  help  to  elliminate  all  but  one  combination. 

A 

> 


j 


A 
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i  . 

VII •  Calculations  in  the  Cubic  Phas'e 

•  *  •  i  '  i 

For  the  cubic’  PbTiO-  crystal  in  the  absence  of  external  fields, 

'  ^  ’  J 

the; only  interesting  observations  are’  the  behavior  of  the  refractive.  > 
index  and  dielectric  ‘constant1 with  temperature.  1 

The  calculated  refractive  index  at  different  temperatures  is  shown 

•  1  i 

in  Figure  7  .  In  this  figure,  the  refractive  index  at  the  curie  temper-* 

i  • 

ature  is  slightly  smaller  than  the  2.72  value  initially  calculated.  . 
This  is  due  to  having  reduced  slightly  the  electronic;  polarizabilities 

i  i 

i  « 

l 

in  an  attempt i to  reduce  the  quadratic  electro^optic  discussed  in  a 

t  1 

later  chapter.  No  experimental  values  of  the  refractive  index  were  < 

•  *  * 

1 

found  for  comparison  at  temperatures  other  than 'the  curie  temperature. 

i 

!  ; 

•  The  dielectric  constant ' s  change  with  temperature  is  'much  larger 

i 

than  the  refractive  index's  change.  •  Calculations  were  done  for  the 

!  i  1 

different  combinations  of  total  polarizabilities  that  satisfied  the. 

i  : 

dielectric  constant  'at  the' cu^ie- temperature.  A  few  of  the  calculated' 

•  .  i 

I  l 

curves  that  best  show  the  effects  of  using  different'  polarizabilities, 

I 

are  plotted  in  Figure  8  .  Since. the  experimental  dielectric 'constant 
versus  temperature  has  been  shown  to  satisfy  the  Curie-Weiss  law 


*  r  J 

for  a  Curie -constant  Q  ='4.11  •  105  '  and  a  Curie-Weiss  temperature 
=  S449°C  (Eef  ■*  6:^-1),  the 'curve- meeting  these  conditions  is  shown 

i 

for  comparison..  '  . 

;  *  1 
Comparison  between  the  theoretical  curves  and  the  Curie-  Weiss  fit  i 

I 

show  that;'  If  the1 theoretical  curve  (#3)  :is  picked  to  fit  at  the  curie 

1  <  ’ 

temperature, •  then  it  falls  far  below  at  higher  temperatures.  And,  if 
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Pig.  7.  Refractive  Index  vs.  Temperature  for  Cubic  P'oTiO^ 
(Theoretical) 
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ihe  theoretical  curve  (#2)  is  raised  towards  the  Curl e-Weiss  curve  at 
higher  tespera tares,  then  it  has  gone  to  infinity  at  the  curie  tenper- 
ature.  Curve  #3  was  the  aost  satisfactory  and,  as  will  he  seen  in  a 
later  chapter,  the  corresponding  polarizabilities  give  the  best  results 
in  the  quadratic  electro-optic  calculations. 
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VIII.  Calculations  in  the  Tetragonal  Phase 

Dielectric  Constants 

!  As  was  discussed  in  Chapter  VI  ,  various  combinations  of  total  polar¬ 
izabilities  were  found  to  produce  the  required  dielectric  constant  along 
the  *a'  axis  at  room  temperature.  Matching  both  the  'a*  and  *c*  axes 
experimental  dielectric  constant,  however,  could  only  be  done  for  a 
combination  that  showed  mixed  positive  and  negative  local  field  factors. 

A  number  of  these  combinations  have  been  used  to  calculate  the  dielectric 
constant  behavior  with  temperature  in  the  tetragonal  phase  to  help  in 
selecting  a  unique  set  for  this  phase.  The  results  of  two  calculations 
are  presented  here. 

The  dielectric  constants,  for  the  polarizabilities  that  matched 
both  axes  at  room  temperature,  along  both  axes  at  various  temperatures 
are  given  in  Figure  9  .  The  *c*  axis  dielectric  constant  shows  a 
gradual  rise  at  low  temperatures,  then  a  more  rapid  rise  and  finally 
rises  very  steeply  at  temperatures  approaching  the  transition  temperature. 
This  is  the  desired  behavior  for  both  axes  (Ref  7).  The  'a'  axis  .s 
unsatisfactory’-,  dropping  throughout  the  temperature  range  rath  no  increase 
even  at  temperatures  approaching  the  transition  temperature.  The  'a' 
axis  dielectric  constant  behavior  wasn't  too  surprising,  as  the  effect 
of  the  mixed  negative  and  positive  local  field  factors  with  temperature 
was  difficult  to  predict. 

The  second  set  of  results  presented  here  are  for  polarizabilities 
with  the  same  titanium  value  as  above,  but  for  lower  lead  and  oxygen  values. 
At  these  lower  values,  the  local  field  factors  are  all  positive  and  the 
results  for  the  'a'  axis  now  show  a  satisfactory  behavior  while  the  'c' 
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Fig.  10.  Dielectric  Constants  vs.  Temperature  in  Tetragonal 
PbTiO^  -  Low  Load  Polarizability  (Theoretical) 
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results  resained  very  low  throughout  (Figure  10  )•  The  low  value  of 
the  'c*  #xis  dielectric  constant  was  characteristic  of  the  calculations  at 
other  polarizabilities.  It  appears  to  be  due  to  a  very  large  effect  of 
the. local  fields  on  the  polarizabilities. 

Covalent-ionic  Bonding  Factors 

Is  was  noted  above,  a  very  large  anisotropy  was  exhibited  by  the  dielec¬ 
tric  constant  in  the  tetragonal  phase  calculations.  This  large  anisotropy 
which  is  partially  produced  by  the  crystal's  internal  electric  fields  via 
the  field-dependent  electronic  polarizabilities,  is  also  exhibited  by  the 
refractive  index.  At  room  temperature,  the  natural  birefringence  was 
calculated  to  be  approximately  0.70  conpared  to  the  experimental  value 
of  approximately  0.01  .  In  this  section,  covalent-ionic  bonding  factors 
are  established  for  reducing  the  local  fields  at  the  ion  sites,  in  an 
effort  to  reduce  the  large  anisotropies. 

The  covalent-ionic  bonding  factor  was  varied  in  two  ways.  The 
bonding  factor  was  first  given  the  same  value  at  all  bonds  and  varied 
from  0.0  to  1 .0  .  Then,  the  factor  was  given  different  values  at  the 
different  bonds.  The  reason  for  treating  the  bonds  individually  was 
because  Wunsch  found  that  two  different  values  of  this  factor  were 
required  to  match  the  experimental  values  of  natural  birefringence  and 
spontaneous  polarization  for  BaTiO^  .  It  was  hoped  that  by  allowing 
different  values  of  covalent  bonding  at  the  different  bonds  (seemingly 
a  truer  physical  situation),  the  correct  birefringence  and  spontaneous 
polarization  might  be  calculated  for  only  one  set  of  covalent  bonding 
factors. 

To  understand  hew  the  bonding  factors  were  determined,  refer  to 
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Figure  11,  of  a  unit  cell  in  the  perovskite  crystal.  If  within  that  unit 
cell,  a  scalier  cell  is  considered,  such  as  that  outlined  by  the  heavy 
lines  with  cne  fall  ion  at  each  lattice  point,  then  the  proper  number 
and  types  of  ions  and  bonds  have  been  accounted  for.  The  different 
bonding  factors  we^-e  determined  by  considering  Shuvaev's  report  (fief  23: 
333)  that  in  BaTiO^  the  minimum  charge  of  the  titanium  ion  is  1.6  units 
of  electronic  charge  equivalent.  With  this  in  mind,  a  charge  of  1.0 
was  assigned  to  the  titanium  in  PbTiOj  based  on  the  comment  that  the 
hi^ier  polarizable  lead  ion  increases  the  bond  covalency  between  titanium 
and  oxygen  (Ref  1:1433).  The  next  step  was  to  give  one-half  of  the 
charge  lost  by  the  titanium  to  the  oxygen  III  ion,  because  of  its  nearness 
in  the  tetragonal  phase,  and  to  share  the  other  half  between  the  oxygen 
I  and  oxygen  II  ions.  The  lead  ion  was  left  with  its  full  charge  of 
+2.0  .  The  resulting  charge  distribution  is  shown  in  Figure  11  .  The 
covalent-ionic  bonding  factors  are  now  simply  equal  to  the  ratio  of  the 
final  charge  of  the  ion  to  its  original  full  ionic  charge. 

A  second  set  of  bonding  factors  was  determined  by  assuming  that  in 
addition  to  the  charge  distribution  described  above,  the  lead  ion  now 
loses  a  fourth  of  its  charge  equally  to  the  oxygen  I  and  oxygen  II 
ions.  The  two  sets  of  bonding  factors  are  listed  in  Table  V  . 
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Pig.  11.  ‘Effective*  Charges  for  Ions  in  Tetragonal 

PbTiO_  (Arbitrary) 

5 


Table  V. 

Special  Covalent-ionic  Bond.'.ig  Sets  for  pbTiO 

_ _ _ 3 

Set  Pb  Ti  Oj  O^ 

1  1.0  0.25  0.625  0.625  0.25 

2  0.75  0.25  0.50  0.50  0.25 
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Natural  Birefringence 

As  already  mentioned  in  the  previous  section,  the  natural  birefrin¬ 
gence  for  the  fully  ionic  crystal  at  room  temperature  was  calculated  to 
be  approximately  0.70  ;  a  value  almost  70  times  larger  than  the  experi¬ 
mental  value  of  0.01  from  Figure  4  .  In  this  section,  the  covalent- 
ionic  bonding  factor  is  employed  in  reducing  the  birefringence. 

At  room  temperature,  the  natural  birefringence  was  calculated  for 
the  various  bonding  factors  discussed  in  the  previous  section:  a  simul¬ 
taneous  variation  from  1.0  to  0.0  at  all  bonds  and  then  the  two 
sets  from  Table  V  .  The  results  are  shown  in  Figure  12  .  The  birefrin¬ 
gence  i'  seen  to  vary  quite  predictably  from  0.198  to  0.68  for 
bonding  factor  values  from  0.0  to  1.0  ,  respectively.  The  zero  bonding 
factor  may  be  physically  unrealizable,  but  it  does  establish  a  minimum 
birefringence  that  may  be  calculated  by  the  use  of  covalent-ionic  bonding 
factors  alone.  The  two  X’s  in  the  figure  correspond  to  the  two  special 
sets  of  factors  from  Table  V  . 

The  natural  birefringence  was  also  calculated  as  a  function  of 
temperature  to  determine  if  perhaps  the  qualitative  behavior  might  be 
correct.  The  results  (Figure  13)  do  not  compare,  as  the  birefringence 
continues  to  increase  after  400°C  ,  while  the  experimental  curve 

(Figure  4)  decreases  above  this  temperature.  This  calculation  corresponds 
to  the  ful'y  ionic  crystal. 


Natural  Birefringence 
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1 - 1 - 1 - 1 - r 


(1) 

(2) 


(Notei  (l)  &  (2)  are  values  at  the  Special  Bonding 
Factor  sets  #1  &  #2  from  Table  V.) 

_ J _ I _ _ I _ I _ 1 _ 


0  0.2  0.4  0.6  0.8  1.0 

Covalent-ionic  Bonding  Factor 

Fig.  12.  Natural  Birefringence  vs.  Covalent-ionic  Bonding 
Factor  for  Room  Temperature  PbTiO^  (Theoretical) 
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Spontaneous  Polarization 

The  spontaneous  polarization  was  calculated  as  a  function  of  the 
covalent-ionic  bonding  factor  at  room  temperature.  In  these  calculations 

!  A 

the  experimental  75/tC/cin  value  (Ref  8ii230)  was  approximately  matched 
by  a  covalent-ionic  bonding  factor  equal  to  0.7  and  by  special  set  #2 
from  Table  V  . 

The  spontaneous  polarization  was  also  calculated  as  a  function  of 
temperature  for  the  full  ionic  crystal  and  is  seen  to  increase  from 
q3/4c/cm^  at  room  temperature  to  88yLfC/cm^  at  470 #C  ,  as  compared  to 
a  reported  value  of  17/4 C/cm^  at  the  transition  temperature  (Ref  I6i44). 

No  attempt  has  been  made  here  to  consider  a  possible  decrease  in  the  ratios 
of  the  ion  shift  distances  to  the  'c'  axis  dimension  with  increasing 
temperature,  as  exhibited  by  BaTiO^  (Ref  19i163)»  since  these  measurements 
have  not  been  reported  for  PbTiO^  .  Doing  so,  plus  treating  the  crystal 
as  only  partialy  ionic,  would  result  in  a  value  closer  to  the  reported 
value  but  would  introduce  undesired  speculation  at  this  point. 

Discussion 

Efforts  to  reduce  the  large  natural  birefringence  and  dielectric 
constant  anisotropy  via  the  covalent-ionic  bonding  factor  have  failed. 

An  interesting  result  is  that  even  for  a  covalent-ionic  bonding  factor 
equal  to  zero  (removing  all  anisotropy  in  the  polarizabilities),  the 
natural  birefringence  is  still  approximately  0.20  .  Since  under  these 
conditions,  onl-y  the  physical  structure,  through  the  lorentz  factors  is 
present  to  produce  these  results,  It  seems  that  this  anisotropy  can  not 
be  reduced  further.  In  Chapter  XI  ,  this  problem  is  considered  by  calcu¬ 
lating  the  contribution  from  each  individual  mechanism  in  the  model. 
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IX.  Electro-optic  Calculations  in  the  Cubic  Phase 
Electro-optic  (g^  -  g^)  Constant 

!  Although  the  model  has  not  been  successful  in  calculating  the  values 
of  the  various  parameters  along  the  *c*  axis  of  tetragonal  PbTiO^  >  it 
has  been  satisfactory  in  the  cubic  phase.  The  model  is  therefore  used 
to  calculate  the  anisotropy  induced  in  the  refractive  index  and  dielectric 
constant  by  a  10  v/M  electric  field.  The  electro-optic  constant  dis¬ 
cussed  here  is  the  quadratic  (g^  -  g^)  constant.  This  constant  has 
been  reported  to  have  a  value  of  approximately  0.01  for  lead-containing 
perovskites  (Ref  24:4790). 

The  elec cro-optic  calculations  were  made  for  the  different  polari¬ 
zabilities  that  predicted  correct  dielectric  constants  at  the  curie 
temperature.  The  results  of  these  calculations  are  presented  in  Figure 
15  .  A  lead  polarizability  approximately  equal  to  that  of  barium  results 
in  the  lowest  value  of  the  quadratic  electro-optic  constant  (curve  #5). 
Going  to  lower  or  higher  values  of  lead  polarizability  increased  the 
constant.  The  lowest  value  (0. 1 3 )  ,  however,  is  an  order  of  magnitude 
larger  than  the  reported  value,  or  approximately  equal  to  the  values  for 
other  barium  titanate-type  perovskites  (Ref  24:4791)  .  The  high  lead 
polarizability  that  had  correctly  predicted  the  • c  *  axis  dielectric 
constant  at  room  temperature,  resulted  in  the  highest  value  for  the 
electro-optic  constant. 

Other  attempts  were  made  to  lower  the  calculated  electro-optic 
constant.  These  attempts  involved:  l)  varying  the  oxygen  total  polari¬ 
zability,  and  2)  varying  the  electronic  polarizabilities  of  titanium 
and  oxygen.  No  indications  were  seen  that  the  0.01  value  could  be 
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approached  by  this  variations. 

There  is  a  strong  possibility  that  the  reported  0.01  value  may 
be  incorrect,  because  it  has  been  calculated  by  using  the  tetragonal 

f 

phase's  natural  birefringence  and  spontaneous  polarization  in  the 
quadratic  electro-optic  relation.  This  assumes  that  such  a  relation, 
the  so-called  "spontaneous"  Kerr  Effect,  holds  in  the  tetragonal  phase. 


The  natural  birefringence  of  PbTiO^  ,  however,  has  not  been  shown  to  be 
proportional  to  the  square  of  the  spontaneous  polarization  as  in  BaTiO^ 
(Ref  19:94),  a~d  this  method  may  not  be  applicable  to  PbTiO^  .  Indeed, 
Remeika  (Ref  16:41 )  has  reported  on  the  incremental  change  of  the 


spontaneous  polarization  with  respect  to  temperature  and  the  spontaneous 


polarization  curve  that  is  implied,  is  not  like  the  birefringence  curve. 


It  also  seems  improbable  that  cubic  PbTiO^  ,  with  crystal  conditions 

similar  to  those  in  cubic  BaTiO^,  should  have  a  quadratic  electro-optic 

effect  so  drastically  different  from  BaTiO,  and  other  perovskites  of 

j 

this  family. 


Dielectric  Constant  Anisotropy 

The  anisotropy  induced  in  the  dielectric  constant  for  cubic  PbTiO., 

j 

has  been  calculated  for  an  electric  field  strength  of  10^  v/M  .  At 
this  field  strength,  the  dielectric  constant  at  the  curie  temperature, 
along  the  applied  field,  drops  to  a  little  less  than  half  its  unbiased 
value.  Although  no  similar  experimental  measurements  are  available  for 
PbTiO^  ,  this  may  be  compared  with  the  results  for  BaTiO^  (Ref  27:116) 
which  showed  a  similar  drop  for  a  field  strength  of  10^  v/M  . 


« 
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X.  Electro-optic  Calculations  In  the  Teti’agonal  Phase 

Since  the  aodel  has  not  been  successful  In  calculating  the  natural 

birefringence  and  dielectric  constants  in  tetragonal  PbTi0o  ,  its  use 
•  j 

in  calculating  electro-optic  effects  in  this  phase  is  questionable. 

For  this  reason,  although  the  calculations  were  done  for  their  instructive 

value,  no  detailed  discussion  is  given.  It  is  sufficient  to  report,  for 

future  considerations,  a  value  for  (r^_  -  (n  /n  )^r1_)  of  -1.03  •  10~^M/V 

jj  o  e  ij 

at  a  field  strength  of  10  v/M  .  In  this  calculations,  the  total  polari¬ 
zabilities  used  were  4.90  A?  ,  0.840  P?  ,  and  2.837  t?  for  lead, 
titanium,  and  oxygen,  respectively.  The  only  value  reported  in  the 
literature  is  a  semi-emperical  value  equal  to  approximately  3  *  10~^m/v 
for  (r^  -  r^j),  (3ef  2O1II6). 
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XI.  Point-dipole  Model  Study 

The  model’s  failure  to  calculate  correctly  the  properties  for  the 
tetragonal  phase  of  lead  titanate  prompted  a  study  of  the  contributions 
to  natural  birefringence  from  the  different  mechanisms  in  the  model. 

The  study  consisted  of  looking  at  the  individual  contributions  from 
tetragonality,  ion  shifts,  and  electric  field  dependence  and  combinations 
of  the  three.  In  doing  the  calculations,  field  effects  were  removed 
by  simply  neglecting  the  field  dependence  of  the  electronic  polarizability, 
while  the  tetragonality  and  ion  shifts  were  considered  by  setting  one 
or  the  other  equal  to  zero. 

The  results  can  be  best  summarized  by  Figure  16  showing  the  different 
contributions  at  different  temperatures.  Several  statements  can  be  made 
from  the  curves  in  the  Figure.  These  are:  a)  Each  individual  contri¬ 
bution  is  larger  than  the  experimental  value.,  b)  The  individual  contri¬ 
butions  all  add  constructively  at  all  temperatures.,  c)  The  shifts 
and  field  dependence  contribution  both  increase  with  temperature.,  and 
d)  Above  the  reference  base  line,  curve  #1  shows  a  behavior  similar  to 
that  of  the  experimental  curve;  a  behavior  which  is  also  evident  in  the 
other  cu  ves.  The  sharp  increase  in  the  topmost  curve  at  temperatures 
approaching  the  transition  temperature  is  due  to  the  sharp  decrease  in 
the  crystal’s  unit  cell  volume. 

It  is  obvious  from  the  results  that  both  the  tetragonality  and  ion 
shifts  are  so  large  in  PbTiO^  ,  that  they  alone  can  produce  a  birefringence 
much  larger  than  the  experimental  value.  One  more  thing  is  clear:  the 
effects  of  tetragonality  and  ion  shifts  add  constructively  and  do  not 
subtract  from  each  other  as  had  been  hoped  for  in  explaining  the 
anomalously  low  birefringence  of  PbTiO,.  . 
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Fig.  16,  Contribution  to  Natural  Birefringence  by  Different 
Mechanisms  in  Point-Dipole  Model 
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XII.  Results.  Conclusions,  and  Recommendations 


Results 


The  application  of  a  point-dipole  model  to  calculate  ferroelectric 


and  electro-optic  properties  of  PbTiO^  has  resulted  in  the  following: 

a)  Satisfactory  refractive  index  and  dielectric  constant  in  the  cubic 
phase  at  the  curie  temperature 

b)  Satisfactory  'a*  ax:is  dielectric  constant  and  spontaneous  polari¬ 
zation  in  the  tetragonal  phase  at  room  temperature 

c)  Acceptable  calculations  of  the  dielectric  temperature  dependence 
in  the  cubic  phase 


d)  An  electro-optic  (g  .  -  g  _)  constant, in  the  cubic  phase,  approxi- 

1 J 

mately  equal  to  that  of  most  barium  titanate-type  perovskites, 
in  contrast  to  a  reported  semi-empirical  value  that  is  only  one- 
tenth  the  magnitude 

e)  An  electro-optic  (r  -  (n  /n  jt.,)  constant,  in  the  tetragonal 

jj  0  6  \j 

—1 1  6 
phase,  equal  to  -1.03  *  10  m/v  for  an  electric  field  of  10  y/'A 


f)  An  unsatisfactory  value  of  the  natural  birefringence  (large  by  a 
factor  of  70  )  and  of  the  anisotropy  in  the  tetragonal  dielectric 
constant 


g)  A  qualitative  behavior  of  the  temperature  dependent  birefringence 
in  the  tetragonal  phase  that  is  not  unlike  the  experimentally 
measured  behavior 


A  study  was  then  made  of  the  different  mechanisms  in  the  model,  by 
calculating  their  individual  contributions  to  the  natural  birefringence. 
The  results  showed  that  disregarding  the  effects  of  the  crystal’s  fields, 
the  crystal's  physical  structure  alone,  through  the  lorentz  factors, 
produced  a  natural  birefringence  that  was  more  th;n  an  order  of  magnitude 
larger  thar  the  experimental  value. 


Conclusions 

The  above  results  lead  to  the  following  reneral  conclusion:  The 
point-dipole  model,  in  the  form  used  by  VJunsch,  can  be  satisfactorily 
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applied  to  PbTiO^  in  its  cubic  phase  to  calculate  ferroelectrics  and 
electro-optic  properties,  and  it  can  not  be  satisfactorily  applied  to 
to  PbTiO-  in  its  tetragonal  phase.  The  unsatisfactory  behavior  of  the 

3  i 

model  in  the  tetragonal  phase  :;eems  to  be  due  to  the  dominating  effect 
of  the  crystal’s  large  c/a  axis  ratio  and  highly  unsymmetrical  ion 
positions. 

Recommendations 

Since  the  present  form  of  the  field-dependent  electronic  polari¬ 
zability  neglects  effects  in  a  direction  perpendicular  to  the  electric 
field,  further  efforts  should  consider  the  significance  of  including 
such  an  effect.  Furthermore,  since  it  is  unlikely  that  the  perpendicular 
effects  alone,  could  correct  the  large  anisotropy  in  refractive  index 
and  dielectric  cons  ant,  further  efforts  should  also  consider  the 
possibility  of  adding  an  explicit  dependence  of  the  polarizabilities  on 
bond  covalency. 

In  the  present  effort,  a  computer  program  was  written  to  operate 
on  the  Air  Force  Weapons  Laboratory's  computer  system.  This  program 
consists  of  the  various  programs  used  by  Wunsch  in  doing  the  various 
basic  calculations  plus  many  added  capabilities,  such  as  temperature 
dependence  and  plotting  capabilities.  This  program  called  FERROEL 
calculates  all  the  desired  calculations  with  a  minimum  number  of 
inputs  and  should  be  used  in  further  efforts  to  improve  the  point- 
dipole  electro-optic  model. 
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Appendix  A 


Electro-optic  Effects 


The  change  in  a  crystal's  refractive  indices  in  the  presence  of  an 
external  electric  field  can  be  best  explained  by  considering  the  equation 
for  the  crystal's  field-free  index  ellipsoid  in  a  non-principal  coordinate 
system: 


I 


(35) 


i 


\ 

i 


$ 


The  presence  of  an  electric  field  results  in  changing  the  coefficients 
of  the  different  terms.  The  change  in  each  coefficient  can  be  expanded 
in  powers  of  the  electric  field.  In  this  expansion,  the  term  that  is 
proportional  to  the  first  power  of  the  electric  field  defines  the 
coefficients  of  the  linear  ( Pockel)  electro-optic  effect.  This  defining 
equation  is 


where  the  p.j  (Pockel  coeff: "ients)  are  members  of  a  second  order 
tensor.  The  expansion  term  that  is  proportional  to  the  square  of  the 
electric  field  defines  the  coefficients  of  the  quadratic  (Kerr)  electro¬ 
optic  effect.  This  defining  equation  is 


RijkLEkEL 


(37) 


where  the  R,  (Kerr  coefficients)  are  members  of  a  fourth  order  tensor. 

ijkL 
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In  the  tetragonal  phase  of  the  perovs’"te  ferroelectrics,  the  linear 
effect  is  much  larger  than  the  quadratic  effect  and  so  the  quadratic 
effect  is  usually  disregarded.  In  the  cubic  phase,  the  linear  effect 

I 

is  not  allowed  because  the  crystal  lattice  has  inversion  symmetry. 

In  a  crystal's  principal  coordinate  system,  the  fourth,  fifth,  and 
sixth  coefficients  of  equation  (35)  are  zero,  so  that  in  the  presence  of 
an  electric  field,  it  may  be  written  as 


Writing  the  A[n2L  in  terms  of  the  Pockel  coefficients  for  the 
tetragonal  perovskite  crystal  (keeping  only  the  non-zero  terms),  for  an 


in  terms  of  the  Pockel  coefficients  for  the 


electric  field  in  the  z  direction,  equation  (38)  becomes  (Ref  27:119) 


(39) 


where  the  new  coefficients  of  the  different  terns  define  the  new  refrac¬ 
tive  indices  along  the  crystal's  principal  axes. 

Subtracting  the  refractive  indices  along  the  x  and  z  directions 
from  each  other,  results  in  the  following  equation  for  the  linear 
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eleciro-optic  constant 


where  the  prime  denotes  values  in  the  presence  of  thb  electric  field. 
Similarly,  for  the  cubic  phase,  equation  (38)  becomes  (Ref  27il00) 


I 


+  R,2E 


+  R,2E 


I 


(41) 


In  this  case,  subtracting  the  new  refractive  indices  along  the  X  and 
Z  directions  from  each  other  results  in  the  following  equation  for  tne 

I 

quadratic  electro-optic  constant . 


2  An 
n3F2 

1  q3 


where  /\n  is  the  induced  birefringence  ih  the  presence  of  the  electric 
field.  .  • 
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